Non-local correlations in the Haldane phase for an XXZ spin-1 chain: A perspective from infinite 

matrix product state representation 
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String correlations are investigated in an infinite-size XXZ spin-1 chain. By using the infinite matrix product 
state representation, we calculate a long-range string order directly rather than an extrapolated string order in a 
finite-size system. In the Neel phase, the string correlations decay exponentially. In the XY phase (Tomonaga- 
Luttinger liquid phase), the behaviors of the string correlations show a unique two-step decaying to zero within 
a relatively very large lattice distance, which makes a finite-size study difficult to verify the non-existence of 
the string order. Thus, in the Haldane phase, the non-vanishing string correlations in the limit of a very large 
distance allow to characterize the phase boundaries to the XY phase and the Neel phase, which implies that 
the transverse long-range string order is the order parameter for the Haldane phase. In addition, the singular 
behaviors of the von Neumann entropy and the fidelity per lattice site are shown to capture clearly the phase 
transition points that are consistent with the results from the transverse long-range string order. The estimated 
critical points including a Berezinsky-Kosterlitz-Thouless transition from the XY phase to the Haldane phase 
agree well with the previous results: A r2 = for the XY-Haldane phase transition and A c3 = 1.185 for the 
Haldane-Neel phase transition from the density renormalization group. From a finite-entanglement scaling of 
the von Neumann entropy with respect to the truncation dimension, the central charges are found to be c - 1.0 at 
A c2 = and c s 0.5 at A rf = 1.185, respectively, which shows that the XY-Haldane phase transition at A c2 = 
belongs to the Heisenberg universality class, while the Haldane-Neel phase transition at A C 2 = 1.185 belongs 
to the two-dimensional classical Ising universality class. It is also shown that, the long-range order parameters 
and the von Neumann entropy, as well as the fidelity per site approach, can be applied to characterize quantum 
phase transitions as a universal phase transition indicator for one-dimensional lattice many-body systems. 

PACS numbers: 75.10.Pq, 75.40.Cx, 75.40.Mg, 03.67.Mn 



I. INTRODUCTION 

Since Landau introducd the theory of second order phase 
transitions, understanding local order parameters character- 
izing different quantum phases has become one of the main 
paradigms in condensed matter physics 1,2 . Local order pa- 
rameters are then known to detect a spontaneous symmetry 
breaking for quantum phase transitions^ 4 -. In some cases, 
however, this local order parameter approach does not work 
because quantum phase transitions arise from the cooperative 
behavior of a system due to emergence of non-local order 5 . 
Although a non-local order parameter could not be directly 
probed in experiments, it can be useful for understanding the 
underlying physics of quantum phases as well as for mark- 
ing phase boundaries. Thus, understanding non-local orders 
in low-dimensional spin systems at absolute zero temperature 
have been an important subject in quantum phase transitions 
in recent years 6 ^— . 

A prototype example, in particular, is the spin-1 antifer- 
romagnetic Heisenberg chain The groundstate of the 
spin chain is in a distinct phase with a finite energy gap, 
but does not exhibit any local order parameter 6 . These 
properties of the groundstate are comprehensively under- 
stood by a non-vanishing non-local string correlation 15 intro- 
duced by den Nijs and Rommelse 16 , and by understanding 
the Affleck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonians for 
interger-spin chains 1718 . The distinct phase of the ground- 
state is called the Haldane phase 19 . The energy gap is also 
called the Haldane gap, which has been manifested by exper- 
imental evidences found in CsNiCL 20 and the organic crystal 



Ni(C 2 H 8 N 2 ) 2 N0 2 C10A 

Indeed, in order to characterize the Haldane phase, such 
non-local string correlations have been extensively studied 
in various finite-size spin systems such as anisotropic spin- 
1 Heisenberg chains 14 , frustrated antiferromagnetic Heisen- 
berg spin-1 chains 18 , alternating Heisenberg chain o 22 ' 23 , spin 
ladders 24-26 and tubes 1 1 , restricted solid-on-solid modeli 6 ., lat- 
tice boson systems 2 - 7 -, and so on. The density matrix renor- 
malization group (DMRG) 28 , the large-cluster-decomposition 
Monte Carlo method 29 , and exact diagonalization with Lanc- 
zos metho d 30 i 31 have been applied for these studies. In 
a recently developed tensor network (TN) representation, 
i.e, matrix product state (MPS) representation 32 , the DMRG 
method 33 also has been applied to explore a string correlation 
for a finite-size lattice. The non-local string order inferred 
from the string correlation behaviors in such finite-size spin 
systems has been used to characterize the Haldane phase from 
other phase s 14,33 . In fact, it is then believed that string correla- 
tions can characterize the Haldane phase. However, no char- 
acterization of the Haldane phase, to the best of our knowl- 
edge, has been made by directly computing long-range string 
order (LRSO) instead of the extrapolated behaviors of string 
correlations till now because all pervious studies have been 
carried out in finite-size systems. 

Thus, in this study, we will investigate string correlations 
and their extreme values for very large spin lattices, i.e., di- 
rectly computing the string order. To do this, we consider 
the infinite-size spin-1 antiferromagnetic Heisenberg chain 
with anisotropic exchange interaction A. We will employ 
the infinite matrix product state (iMPS) representation^— 



for the ground state wavefunction of the infinite lattice sys- 
tem. The groundstate wavefunction can be obtained numer- 
ically by using the infinite time evolving block decimation 
(iTEBD) method 35 within the iMPS representation. To inves- 
tigate string correlations, we will introduce an efficient way 
to calculate a non-local correlation and its extreme value for 
a large-size lattice system in the iMPS representation. It is 
found that, except for the Haldane phase, the string correla- 
tions decay exponentially in the Neel phase, while they show 
a unique behavior of decaying to zero within very large lattice 
distance in the XY phase (Tomonaga-Luttinger liquid phase). 
For the Haldane phase, the string correlations are saturated 
to finite values, which shows a LRSO as the lattice distance 
goes to infinity. Also, from the LRSO with respect to the 
anisotropic interaction parameter A, it is clearly shown that 
both the x- and y-components rather than the z-component of 
the LRSO play a role as the order parameters characterizing 
the Haldane phase. As a consequence, the string order param- 
eters enable us to directly characterize the possible phases of 
the system with respect to the anisotropic exchange interac- 
tion. Moreover, the von Neumann entropy and the fidelity per 
lattice site (FLS) are calculated to show that their singular be- 
haviors correspond to the phase transition points. The central 
charges from the finite-entanglement scaling quantify the uni- 
versality classes of the transition points. The FLS is shown to 
capture a Berezinsky-Kosterlitz-Thouless (BKT) type transi- 
tion, in contrast to the fidelity susceptibility that fails to detect 
it. 

This paper is organized as follows. In Sec. II, a brief expla- 
nation for the iMPS representation is given. We discuss how 
to capture non-local correlations including string correlations 
and string orders directly by exploiting the iTEBD method. In 
Sec. Ill, the spin-1 XXZ chain model is introduced. We dis- 
cuss the behaviors of the string correlations and Neel correla- 
tions as a function of the lattice distance for given anisotropic 
interaction strengths in Sev. IV. The phase diagram of the 
spin-1 XXZ chain model is presented based on the non-local 
correlations and the string and Neel order parameters in Sec. 
V. In Sec. VI, we discuss local and non-local properties of the 
iMPS groundstate that allow to introduce pseudo symmetry 
breaking order for the XY phase for finite truncation dimen- 
sions. In Sec. VII, the phase transitions and their universality 
classes are discussed from the von Neumann entropy and the 
central charges via the finite-entanglement scaling. In Sec. 
VIII, the groundstate FLS is shown to have a clear pinch point 
that corresponds to a quantum phase transition. In Sec. IX, 
our conclusions and remarks are given. 
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FIG. 1 : (color online) (i) Diagrammatic infinite matrix product state 
(iMPS) representation of a wavefunction I*?) having two-site trans- 
lational invariance for infinite one-dimensional lattice systems. A A ^ 
denoted by filled diamonds are a diagonal (singular value) matrix, 
respectively, depending on the 2i-th (A) and 2i + 1-th (B) bonds. 
r A(B) indicated by filled circles are a three-index tensor for the 2i- 
th and 2i + 1-th sites, respectively, r and / indicate the bond indices, 
(ii) A compact form of the pictorial representation in (i), where A 
and B (filled circles) by absorbing A into the tensors T (indicated by 
the dashed-line box in (i)) present the lattice sites, (iii) The norm 
( V P| V P) of a given state \T) in (i). The black dots indicate the left and 
right dominant eigenvectors that can be determined for the dominant 
eigenvalue of the transfer matrix T. The transfer matrix T is obtained 
by the contraction on the tensors A, A", B, and B* in the dashed-line 
box. 



function represented in TNs allows to perform the classical 
simulation of quantum many-body systems. Especially, in 
one-dimensional spin systems, a wave function for infinite- 
size lattices can be described by the iMPS representation 32 . 
The iMPS representation have been successfully applied to in- 
vestigate the properties of ground-state wave functions in var- 
ious infinite spin lattice systems. The examples include Ising 
model in a transverse magnetic field 36 and with antisymmetric 
anisotropic and alternative bond interactions 37 , XYX model 
in an external magnetic field 36 , and spin-1/2 XXZ model 38 . 
However, the iMPS has not been applied much to explore spin 
correlations. Few studies have shown the behaviors of spin- 
spin correlations in the infinite Ising spin chain 35 . Further- 
more, non-local spin correlations have not been explored yet 
in infinite-size systems with the iMPS representation. Then, 
in this section, we will discuss how to calculate a non-local 
spin correlation within the iMPS representation. 



A. iMPS representation and iTEBD algorithm 



For an infinite one-dimensional lattice system, a state can 
be written as 35 ^ 



II. IMPS REPRESENTATION AND NON-LOCAL 
CORRELATIONS IN NUMERICAL METHOD 

Recently, significant progress has been made in numerical 
studies based on TN representations 3 ^^ for the investigation 
of quantum phase transitions, which offers a new perspec- 
tive from quantum entanglement and fidelity, thus providing 
a deeper understanding on characterizing critical phenomena 
in finite and infinite spin lattice systems. Actually, a wave 
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where \S i) denote a basis of the local Hilbert space at the site 
i, the elements of a diagonal matrix are the Schmidt de- 
composition coefficients of the bipartition between the semi- 
infinite chains L(-oo A and R(i+ 1, oo), and F [/ ' e are 
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a three-index tensor. The physical indices St take the value 
1, • • ■ , d with the local Hilbert space dimension d at the site i. 
The bond indices a, take the value 1, ■ • • ,x w i m the trunca- 
tion dimension of the local Hilbert space at the site i. The 
bond indices connect the tensors F in the nearest neighbor 
sites. Such a representation in Eq. (Q~|) is called the iMPS 
representation 35 . If a system Hamiltonian has a translational 
invariance, one can introduce a translational invariant iMPS 
representation for a state. Practically, for instance, for a two- 
site translational invariance, the state can be reexpressed in 
terms of only the three-index tensors T A{B) and the two diag- 
onal matrices Aa(b) for the even (odd) sites 37 , where {T, A] are 
in the canonical form, i.e., 
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v, here I and r are the left and right bond indices, respectively. 
In Fig. [T](i), a state I^F) with a two-site translational invariance 
is pictorially displayed in the iMPS representation for infinite 
one-dimensional lattice systems. In a more compact form, fur- 
ther, the quantum state can be reexpressed as the state in Fig. 
[T](ii) by absorbing the diagonal matrices A into the tensors F. 

Once a random initial state |*F(0)) is prepared in the iMPS 
representation, one may employ the iTEBD algorithm^ 5 - to 
calculate a groundstate wavefunction numerically. For in- 
stance, if a system Hamiltonian is translational invariant 
and the interaction between spins consists of the nearest- 
neighbor interactions, i.e., the Hamiltonian can be expressed 
by H - J^fh^' l+1 \ where h! f ' l+1 ^ is the nearest-neighbor two- 
body Hamiltonian density, a groundstate wavefunction of the 
system can be expressed in the form in Eq. (ffj). The imaginary 
time evolution of the prepared initial state |*F(0)), i.e., 
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leads to a groundstate of the system for a large enough r. By 
using the Suzuki-Trotter decomposition 43 , actually, the imagi- 
nary time evolution operator U = exp[-i/r] can be reduced to 
a product of two-site evolution operators U(j, i + 1) that only 
acts on two successive sites i and i+ 1 . For the numerical imag- 
inary time evolution operation, the continuous time evolution 
can be approximately realized by a sequence of the time slice 
evolution gates U(i, i + 1) = exp [-/i [m+1] c)t| for the imagi- 
nary time slice 6t = t/« <k 1. A time-slice evolution gate 
operation contracts Fa, Tb, one Aa, two Ab, and the evolution 
operator U(i, i + 1) = exp [-/! [, ' /+1] c)t|. In order to recover 
the evolved state in the iMPS representation, a singular value 
decomposition (S VD) is performed and the x largest singular 
values are obtained. From the SVD, the new tensors Fa, F#, 
and Aa are generated. The latter is used to update the tensors 
Aa as the new one for all other sites. Similar contraction on 
the new tensors F^, Tb, two new Aa, one Ab, and the evolution 
operator U(i + 1, i + 2) = exp [-/i [,+u+2] ], and its SVD pro- 
duce the updated F^, Tb, and Ab for all other sites. After the 
time-slice evolution, then, all the tensors F^, T B , A A , and A B 
are updated. This procedure is repeatedly performed until the 
system energy converges to a groundstate energy that yields 
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FIG. 2: (color online) Diagrammatic representations of the expec- 
tation values in the iMPS representation. The black dots in left 
(right) side denote the dominant eigenvectors of the transfer ma- 
trix. In between a given wave function IT) and its complex conju- 
gate, each spin operator act on each site, (i) Local order compu- 
tations including one, two, and three spins such as magnetization, 
dimer order, and chiral orders. Non-local correlations can also be 
calculated in the iMPS representation, e.g., (ii) string correlation 
(r s (i,j) = -<57exp[OT2£! +1 5, ff ]5;>(« 6 {x,y,z\\ (iii) parity cor- 
relation Op(i, j) = (exp[(7r2j = , 5^]), and (iv) Neel correlation (two 
point spin correlation) O^(iJ) '= (-l) H (S^S i p. 



a groundstate wavefunction in the iMPS representation. The 
normalization of the groundstate wavefunction is guaranteed 
by requiring the norm < V P| V P> = 1 in Fig. [T](iii). 



B. Non-local correlations 

In principle, once one obtains a groundstate wavefunction, 
the expectation values of physical quantities can be calculated. 
In Fig. [2] we depict the diagrammatic iMPS representations 
for some examples of various expectation value calculations. 
Figure [2] (i) presents the computation of successive spin op- 
erators such as magnetization 36 (Sf), dimer order 44 

and chiral orderi 5 ^ (SfS p i+1 S 7 i+2 ) (a,0,y e {x,y,z}). On cal- 
culating the expectation values, each spin operator acting on 
a site is sandwiched between a given wave function I 1 ?) and 
its complex conjugate. The left and right dominant eigenvec- 
tors of the transfer matrix denoted by the black dots act on 
the tensor, for instance, contracting the tensors A, A*, (or B, 
B*) and the local spin operator Si for (S"). This leads to the 
expectation value (S"). 
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FIG. 3: (color online) String correlation j) for the Haldane spin- 
1 chain with truncation dimension ;^ = 32 in the iMPS representation. 
All components of the string correlations are well saturated to the 
value O a s = 0.37434447, which agrees very well with the result from 
the DMRG in Ref. El 



\i - j\ with the truncation dimension x — 32. It is shown 
clearly that the 0£{i,j) starts to saturate around the lattice 
distance \i - j\ a 10 and Og(i,f) = 0*(i,f) = for 
the Haldane chain. The saturated value of the string cor- 
relation from our iMPS representation approach is given as 
(9g(oo) = 0.37434447, which agrees very well with the value 
Cf^oo) - 0.37432509 as well as the saturation behavior of 
the string correlations from the DMRG method in Fig. 5 of 
Ref. H^. As is well-known, further, the spin-1 AKLT model 
H^At = Z/[(l/2)Si • S /+ i + (1/6)(S, • S, + l ) 2 + 1/3] is ex- 
actly solvable and the string order is given as the exact value 
4/ 9 - 17 ' 18 . In our iMPS representation, the value of the string or- 
der has been confirmed to be 4/9 for the spin-1 AKLT model 
within the machine accuracy. 



III. SPIN- 1 XXZ HEISENBERG CHAIN 



This procedure can be simply expanded to the calculation 
of spin-spin correlations as well as non-local correlations. Ex- 
amples are the string correlation Os 16 , the parity correlation 
Op 46 , and the Neel correlation O^ 44 given by, respectively, 
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where i and j denote the site locations in the lattice and then 
the lattice distance is \i-j\. Figures0(ii), (iii), and (iv) present 
the computation of the string, parity, and Neel correlations, re- 
spectively, in the iMPS representation. Note that, for the cal- 
culation of these correlations, the (z - j) tensors are involved. 
For the parity correlation, the multi-site operator e\p[in £ S£] 
acts on each site between the site i and site j. For the string 
correlations, the spin operators 5,- and S j act on sites i and 
j while the multi-site operator explin^S^] acts on the sites 
in between the sites i and j. Compared to the string corre- 
lations, the Neel correlations can be calculated by replacing 
the multi-site operator expfzVr £ S?] with the identity operator, 
i.e., 7, + i ■ • • 7j_i while the parity correlation can be obtained by 
expanding the multi-site spin operator to the both ends of the 
lattice distance. Then, it should be noted that, in principle, 
the iMPS representation allows to calculate any correlations 
in the limit of the infinite distance, i.e., \i - j\ — > 00. For nu- 
merical calculations, in order to obtain the correlations in the 
limit of the infinite distance, one can set a truncation error s 
rather than the lattice distance, i.e., O a (i, j) - O a (i, j + 1) < s. 
In this study, for instance, s = 10~ 8 is chosen. 

As an example, for the spin-1 Haldane chain H = 2; S, ■ 
S, + i, the exact-diagonalization calculations of a 14-site lat- 
tice have estimated 0|(oo) =2 0.38^. The DMRG methods 
have verified the existence of LRSO by estimating O^(oo) = 
0.37432509 for the Haldane chain* 8 .. In Fig. [3] we plot the 
string correlation (i, j) as a function of the lattice distance 



Spin-1 Heisenberg chains are one of the prototypical ex- 
amples in understanding non-local correlations 8 , i.e., string 
correlation. Then, to investigate non-local correlations in one- 
dimensional spin systems, we consider an infinite spin-1 XXZ 
Heisenberg chain described by the Hamiltonian 

OO 

H = jY,[S*S* +l+ S*Sl 1+ AS*Sl 1 l (5) 



where Sfia = x,y,z) are the spin-1 operators at the lat- 
tice site i, J denotes the antiferromagnetic spin-exchange 
interaction between the nearest neighbor spins, and A is 
responsible for the anisotropy of the exchange interaction. 
This model has been intensively studied for a couple of 
decades ^ 8 ' 28-3 1 i 33 i 49 ^ . The studies have shown that there are 
the four characteristic phases with respect to the anisotropic 
exchange interaction A. If the anisotropic interaction is much 
smallerthan -1, i.e., A <K -l,the Hamiltonian can be reduced 
to a spin-1 ferromagnetic Ising model H w - Tji S^S z l+l and 
the system is in the ferromagnetic phase. At A c i = -1, a first- 
order transition occurs between the ferromagnetic phase and 
the XY phase. If the anisotropic interaction is much greater 
than 1, i.e., A » 1, the Hamiltonian is reduced to a spin-1 an- 
tiferromagnetic Ising model H as YiiS^S 1 -,^ and the system is 
in the antiferromagnetic (AF) phase. AtA C 3 = 1.17 + 0.02, the 
Haldane-Neel phase transition occurs, which belongs to the 
two-dimensional Ising universality class^i22. In between the 
two phase transition points A c j < A < A C 3, the XY-Haldane 
phase transition occurs, which has been thought to be a BKT 
transition at A c2 ^— . 

As is well-known, many antiferromagnetic spin systems 
can be explored by the standard spin-spin correlations (Neel 
correlations) 44 . However, in the Haldane phase, the spin- 
spin correlations decay exponentially with a finite correlation 
length 3348 ' 49,58 and the Haldane gap exists. In this aspect, 
the Haldane phase can be considered as a disordered phase 12 . 
Also, the XY phase is characterized by the power law decay 
of the spin-spin correlations (Neel correlations) 58 with gapless 
excitations. Characterizing both the XY and Haldane phases 
therefore is a non-trivial task in the aspect of the spin-spin 
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correlations. By investigating the spin correlations, the transi- 
tion point A C 2 has been estimated to be < A C 2 < 0.2 from 
the exact numerical calculations and the finite-cell-scaling 
analysis 30 , A C 2 = -0.01 + 0.03 from the phenomenological 
renormalization-group technique and the finite-size scaling 
analysis with 16 spin sites 28 , and A C 2 = 0.068 + 0.003 from 
the criterion exponents of the spin correlations rj x = 1 /4 in the 
exact diagonalization method with 16 spin sites 31 . Investigat- 
ing the excitation gap, as the anisotropic interaction strength 
varies, is also a method to characterize the Haldane phase. 
By using the lowest-energy levels and finite-size scaling for 
energy gaps from the Lanczos method, the critical point has 
been conjectured to be at A C 2 = in Ref. |8j As an alternative 
way to characterize the Haldane phase, the string order aris- 
ing due to the fully broken Z2 x Z2 hidden symmetry has been 
investigated 67 . The XY-Haldane transition point has been es- 
timated A C 2 ~ by exploring the string correlations from a 
finite size analysis 14 and by a finite-size scaling of the string 
order—. 

IV. BEHAVIORS OF THE STRING AND NEEL 
CORRELATIONS IN SPIN-1 XXZ CHAIN 

A non-vanishing correlation in the limit of the infinite lat- 
tice distance (\i - j\ — » 00), i.e., a long-rang order reveals that 
the system is in a ordered state. For the non-local correlations, 
the string and Neel orders are respectively defined by 

0? = lim Og(iJ), (6a) 

|!-j|-»00 

O* = lim 0°{i,j). (6b) 
|i-jl-»<» 

For instance, the non-vanishing spin-spin (Neel) correlations 
for \i - j\ — > 00 indicate that the system is in an antiferro- 
magnetic state. Also, the ground state in the Haldane phase is 
known to be characterized by the string order 7 . In the view- 
point of the string order, then, the Haldane phase could be an 
ordered phase. Further, if the string order plays a role as the 
order parameter for the Haldane phase, from the string order, 
the phase transition boundary from the Haldane phase to other 
phases can be captured. This view has been applied to the in- 
vestigations of the Haldane phase of spin-1 systems^ - - — . 
By using numerical-diagonalization, however, available sys- 
tem sizes were too small to convince string order behaviors as 
an order parameter. Thus, in Ref. |33l comparisons between 
the behaviors of the string and the Neel correlations from a 
finite size spin lattice (up to 300 sites), and their finite-size 
scaling behaviors have been used to capture the phase bound- 
ary. However, directly capturing the critical behavior of the 
string order near the transition point was quite difficult due to 
a very limited lattice size. Compared to such approaches, as 
discussed in Sec. II, the iMPS approach enables to explore the 
behaviors of the string order directly in the limit of the infinite 
lattice distance (|z - j\ — » 00). 

In Fig. |U we plot the string correlations 0$ (i, j) as a func- 
tion of |/ - j'| for various anisotropic interactions A. In Figs. 
(c) for the Haldane phase and (d) for the Neel (antiferro- 
magnetic) phase, the string correlations show a logarithmical 
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FIG. 4: (color online) String correlations as a function of the 
lattice distance with the truncation dimension ;^ = 32 for various 
anisotropic interaction (a) A = -0.2 and (b) A = 0.364 in the XY 
phase, (c) A = 0.6 in the Haldane phase, and (d) A = 1.3 in the 
Neel phase. In fact, for the truncation dimension x = 32, the phase 
transition points are given as A C 2Ct = 32) = 0.366 and A c i(x = 32) = 
1.180 in Fig. [6] It should be noted that, in the XY phase, the string 
correlations decay to zero with a unique two-step decaying behavior. 



decaying to its saturated value or zero as the lattice distance 
|; - j\ increases up to a few hundreds. While, in the XY phase 
in Figs. 0(a) and (b) , the string correlations show a unique 
two-step decaying to zero. As the lattice distance \i - j\ in- 
creases, that is, the string correlations undergo a decaying be- 
havior for a few hundreds of the lattice distance, a saturation- 
like behavior for a few thousands of the lattice distance, and 
then eventually decaying again down to zero around a few 
tens of thousands. Hence, in contrast to the Haldane phase, 
there is no long-range string order in the XY phase. Also, 
it should be noted that, near the transition point in the XY 
phase in Fig. |4](b), such saturation-like behaviors of the string 
correlations occur for a very wide range of the lattice dis- 
tance from a few hundreds to a few thousands, i.e, roughly 
10 2 \i — j\ < 10 4 . While, away from the transition point in 
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FIG. 5: (color online) Neel correlations 0^(i, j) in the Haldane phase 
(A = 0.8) with truncation dimension ^ = 32. All components of the 
Neel correlations decay to zero exponentially. The transverse (x- and 
y- components) Neel correlations have a same value, i.e., = O y N . 



Fig. 0(a), such saturation-like behaviors of the string correla- 
tions occur for a relatively narrow range of the lattice distance 
roughly 2 x 10 2 < \i — j\ < 10 3 . This saturation-like behav- 
ior of the string correlations makes the characterization of the 
Haldane phase quite difficult directly from a string order be- 
havior in finite size systems. Indeed, compared to our iMPS 
results, a finite size system in the XY phase has given a finite 
value of the string correlation 0^(i, j) within the limitation of 
the finite system sizes (up to 300 sites) 33 . 

In Fig. [5] a Neel correlation is displayed as a function of 
the lattice distance |/ - j\. It is shown that, in the Haldane 
phase A = 0.8, the spin correlations (Neel correlations) de- 
cay exponentially to zero, which allows to characterize the 
phase transition from the Neel phase to the Haldane phase. 
As shown in Fig. 0(d), in the Neel phase, in contrast to the 
z-component of the string correlations that survives for very 
large distances, the x- and y-components of the string corre- 
lations also decay exponentially to zero. Also, as shown in 
Fig. 0(c), all the components of the string correlations in the 
Haldane phase have non-zero values in the limit of the infinite 
lattice distance. Then, alternatively, the x- and y-components 
of the string order make it possible to distinguish the Haldane 
phase from the Neel phase. Further, in the XY phase in Fig. 
0(b), all the components of the string order become zero with 
the unique two-step decaying behavior. As a consequence, the 
x- and y -components of the string order play a role of a true 
order parameter characterizing the Haldane phase from the 
XY phase and the Neel phase. 



V. HALDANE PHASE AND ORDER PARAMETER 

As discussed in Sec. IV, the LRSOs can characterize the 
Haldane phase. In Fig. [6] we plot (a) the string orders 0% 
and (b) the Neel order parameter 0* as a function of the 
isotropic exchange interaction strength A for the truncation 



dimension x = 32. It is shown that the Ol 



have 



non-zero values for 0.366 < A < 1.180, while the <9j has 
a finite value for A > 1.180. The string orders become 
zero for -1 < A < 0.366. Also, the Neel order parame- 
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FIG. 6: (color online) (a) String order parameters 0^' y for the Hal- 
dane phase and (b) Neel order parameter for the Neel phase as a 
function of A with truncation dimension^ = 32. Note that the string 
order O z s does not play a role as the order parameter to characterize 
the Haldane phase in (a). For this truncation dimension^ = 32, the 
system is in the XY phase for -1 < A < 0.366, the Haldane phase for 
0.366 < A < 1.180, and the Neel phase for A > 1.180. For the trun- 
cation dimension^ = 32, then, the transition points are A r 2 = 0.366 
and A c3 = 1.180. 



ter : N has a non-zero value for A > 1.180, which charac- 
terizes the Neel phase. This implies that the Oj and O y s are 
the order parameters for the Haldane phase. Then, the Hal- 
dane phase exists in the range of the anisotropic interaction 
strength 0.366 < A < 1.180, which implies A C 2 = 0.366 and 
A f 3 = 1.180 for the truncation dimension^ = 32. Thus, the 
XY phase occurs for -1 < A < 0.366. 

Actually, the transition points between the phases depend 
on the truncation dimension^-, i.e., A C 2 = A^CtO and A C 3 = 
A c i(x)- As the truncation dimension^ increases from a lower 
truncation dimension (e.g., x — 8), the phase transitions 
A c i(x) an d A c i(x) occur starting at the lower and the higher 
values of A's, respectively. Then, the critical points A C 2(oo) 
and A e 3(oo) in the thermodynamic limit can be extrapolated to 
X — > oo. In Fig. |7j we plot the transition points (a) A C 2(x) 
and (b) A c j,(x) as a function of the truncation dimension^. 
We employ an extrapolation function A(x) — a + bx~ c , char- 
acterized by the fitting constants a, b, and c, which guaran- 
tees that A(oo) becomes a finite value. The numerical fittings 
give, respectively, a = 0.015, b = 1.366 and c = 0.387 for 
the phase transition between the XY and Haldane phases and 
a = 1.185, b = -1.748, and c = 1.800 for the phase tran- 
sition from the Haldane phase to the Neel phase. In Fig. [7] 
in the limit of the infinite truncation dimension, i.e., x — > °°. 
the fitting function are shown to saturate well to the extrap- 
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FIG. 7: (color online) Phase transition points A^Cf) and & c t,(x) as a 
function of the truncation dimension x- Here, the truncation dimen- 
sions are taken as^ = 8, 16, 25, 32, 50, and 64. For the numerical 
extrapolation, the fitting functions are chosen as A c (^) = a + b\~ c , 
where a, b, and c are a real number. A best numerical fitting gives (a) 
a = 0.015, b = 1.366, and c = 0.387 for the XY-Haldane phase tran- 
sition A c2 (x), and (b) a = 1.185, b = -1.748, and c = 1.800 for the 
Haldane-Neel phase transition A c3 (x). The estimated critical points 
in the thermodynamic limit are given by A c2 (oo) = 0.015 for the XY- 
Haldane phase transition and A r 3(oo) = 1.185 for the Haldane-Neel 
phase transition. The critical points agree well with the results from 
the previous studie a 28 ~ 31 t 33 . 



olated value A(oo) = a which can be regarded as a critical 
point A e = A(oo). As a result, our extrapolations give the 
critical points A C 2(°o) = 0.015 and A^oo) = 1.185. Our crit- 
ical points agree well with the results A C 2 = 0.068 ± 0.003 
from the exact diagonalization 31 ., A C 2 — -0.01 + 0.03 from the 
phenomenological renormalization group with the finite-size 
scaling analysis^, A e3 = 1.17 + 0.0222-.2P. and A f3 = 1.18622 
from the DMRG. 
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FIG. 8: (color online) (a) Transverse Neel order 0' N = + O y N as 
a function of the anisotropic interaction A with for various values of 
the truncation dimension x- Note that, in the XY phase, the trans- 
verse Neel order is not zero for the finite truncation dimensions, (b) 
Transverse Neel order 0' N (x) as a function of the truncation dimen- 
sion^ for A = 0.0, -0.5, and -0.7. For the numerical extrapolation, 
the fitting functions are chosen as 0' N (x) = a + bx~ c , where a, b, 
and c are a real number. A best fitting gives the fitting constants as 
(i) a = 0.446 x 10~ 4 , b = 0.948, and c = 0.315 for A = 0.0, (ii) 
a = 3.643 x 10~ 3 , b = 0.936, and c = 0.203 for A = -0.5, and (iii) 
a = 0.897 x 10"*, b = 0.956, and c = 0.156 for A = -0.7. These 
results show that, in the thermodynamic limit, the transverse Neel 
order does not exist in the XY phase, i.e., 0' N (x = oo) = 0, 



VI. XY PHASE 

The XY phase is known to have a power-law decay of the 
spin-spin correlations with a gapless excitation. This implies 
that there exists no long-range order in the XY phase in the 
thermodynamic limit. Actually, in numerical approaches, di- 
rectly characterizing a XY phase from a power-law decay of 
the spin-spin correlation is a non-trivial task. This may be the 
reason why a level spectroscopy of numerical approaches has 
been invented as a useful way to characterize the XY phase. 
However, directly detecting a vanishing excitation gap from 
numerical calculations near the transition point is also not a 
trivial work due to a very limited lattice size. As discussed in 
Sec. V, our iMPS approach has verified that the transverse 
string order parameter and the longitudinal Neel order pa- 
rameter clearly characterize the Haldane phase and the Neel 
phase, respectively, in Figs. 0(a) and (b). Obviously, the lon- 
gitude Neel order and all the components of the string order 
become zero in the XY phase even for the finite truncation 
dimension ;f. 



However, such vanishing behaviors of the long-range order 
parameters in the XY phase do not guarantee the existence of a 
XY phase in the iMPS representation. In this sense, it would 
be worthwhile to discuss a practical way to characterize the 
XY phase in the iMPS representation by defining a pseudo 
order parameter for a finite truncation dimension x- Then, 
some local or non-local properties of our iMPS groundstate 
will be introduced as indicators that can be used to distinguish 
a XY phase from other phases in the iMPS representation22iS. 



A. Transverse Neel order for finite truncation dimensions 

Let us consider the transverse Neel order in the XY phase. 
In Fig. [8] (a), we plot the transverse Neel order as a function 
of the anisotropic interaction strength A for various trunca- 
tion dimensions. Here, the transverse Neel order is defined 
by the sum of the x- and y- components of the Neel order, 
0' N (x) - Oi(x) + O y N (x). It is shown that the transverse Neel 
order has a finite value in the XY phase. Actually, as the 
truncation dimension increases, as discussed in Sec. V, the 
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interaction parameter range of the XY phase becomes nar- 
rower because the transition point A C 2(X) between the XY 
phase and the Haldane phase moves to a lower value for a 
higher truncation dimension. Note that, from the transverse 
Neel order, the transition points between the XY phase and the 
Haldane phase (non-zero values of the transverse Neel order) 
are the same with the values from the string order parameter, 
while the ferromagnetic-XY transition point does not move at 
A c i = -1.0. Thus, the non-vanishing transverse Neel order 
can be used as a pseudo order parameter characterizing the 
XY phase for a finite truncation dimension^-. 

Also, it should be noted that the overall amplitude of the 
transverse Neel order becomes smaller as the truncation di- 
mension increases in Fig. [8] (a). In order to understand 
the transverse Neel order in the thermodynamic limit, i.e., 
X — > oo, in Fig. [8](b), we plot the transverse Neel order as a 
function of the truncation dimension ^ for, as examples, three 
anisotropic interaction strengthes A = 0, -0.5, and -0.7. It 
is shown clearly that the transverse Neel order decreases as 
the truncation dimension increases. We perform an extrapo- 
lation by especially introducing a power-law fitting function 
0' N (x) - a + bx~ c with respect to the truncation dimension^. 
The numerical fittings give (i) a = 0.446 x 10~ 4 , b = 0.948 
and c = 0.315 for A = 0, (ii) a = 3.643 x 10~ 3 , b = 0.936 
and c = 0.203 for A = -0.5, and (iii) a = 0.897 x 10~ 4 , 
b = 0.956 and c = 0.156 for A = -0.7. Hence, the extrapo- 
lated values of the transverse Neel order for ^ — » oo are given 
as 0' N (oo) = 0.446 x 10~ 4 , 3.643 x 10~ 3 , and 0.897 x 10~ 4 
for A = 0, -0.5, and -0.7, respectively. This implies that, 
similar to the power-law decay of spin-spin correlations with 
respect to the lattice distance, the transverse Neel order fol- 
lows a power-law decaying to zero with respect to the trunca- 
tion dimension x- These results show that, in the XY phase, 
the transverse Neel order as well as the longitudinal one also 
becomes zero, 0£(oo) = 0. As a consequence, both the string 
and the Neel long-range orders do not exist in the XY phase. 



B. Pseudo local order for finite truncation dimensions 

Recently, a pseudo local order have been suggested for the 
XY phase in the iMPS representation in Ref. |38l The local or- 
der can be defined as Ol = y(S x ) 2 + (S v ) 2 . In Fig. |9](a), we 
plot the pseudo local order Ol as a function of the anisotropic 
interaction strength A for various truncation dimensions. It is 
shown that the defined local order has a finite value in only 
the XY phase. Similar to the transverse Neel order, from 
the defined local order, the transition points between the XY 
phase and the Haldane phase (non-zero values of the defined 
local order) are detected at the same values from the string 
order parameter for the same truncation dimension x- Also, 
the ferromagnetic-XY transition points from the defined local 
order does not move at A c i = -1 .0. Compared with the trans- 
verse Neel order, Fig. 0(a) shows that the overall amplitude 
of the defined local order becomes smaller as the truncation 
dimension^ increases. In Fig. [9] (b), we plot the transverse 
Neel order as a function of the truncation dimension x for, 
as examples, three anisotropic interaction strengthes A = 0, 
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FIG. 9: (color online) (a) Pseudo local order L = ^(S x ) 2 + (S y ) 2 
for the XY phase as a function of the anisotropic interaction A for 
various values of the truncation dimension x- (b) Pseudo local order 
Ol(x) as a function of the truncation dimension^ for A = 0.0, -0.5, 
and -0.7. By using the fitting function Otix) = a X + C X ') with 
a, b, and c being a real number, the numerical extrapolations give the 
fitting constants as (i) a = 0.980, b = 0.159, and c = -0.028 for 
A = 0.0, (ii) a = 0.998, b = 0.108, and c = -0.120 for A = -0.5, and 
(iii) a = 0.967, b = 0.075, and c = 0.048 for A = -0.7. 



-0.5, and -0.7. It is shown clearly that the defined local order 
decreases as the truncation dimension increases. We perform 
an extrapolation with respect to the truncation dimension^- by 
using the same fitting function Ol(x) = +cy~'), with 

a, b, and c being a real number, given in Ref. [3i| Figure 
[9] (b) shows the behaviors of the defined local order for the 
XY phase agree well with the results of spin- 1/2 XXZ model 
in Ref. (38j Hence, similar to the transverse Neel order, the 
defined local order can be used as a pseudo order parameter 
characterizing the XY phase for a finite truncation dimension 
X- 



VII. ENTANGLEMENT ENTROPY, CENTRAL CHARGE, 
AND UNIVERSALITY CLASS FOR PHASE TRANSITIONS 

Instead of using order parameters, recently, various types 
of quantum entanglement measures have been proposed as an 
indicator characterizing quantum phase transitions 2,60 . One of 
successful measures is the von Neumann entropy for a bipar- 
tite system 61-64 . Singular behaviors of bipartite entanglements 
for a pure state reveal quantum critical behaviors, which has 
been verified as being universal by extensive studies in many 
one-dimensional systems 6 ^ 6 -. 
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FIG. 10: (color online) Von Neumann entropy as a function of the 
anisotropic interaction A for the truncation dimension x = 32. The 
transition points are seen at A c i = -1.0, A c i = 0.366, and A r 3 = 
1.180. This result is consistent with the phase transition points from 
the order parameters in Fig. [6] 



Von Neumann entropy singularities.— In the iMPS ap- 
proach, the von Neuman entropy can be explored. Let us 
recall the diagonal matrix A. As discussed in Sec. II, the el- 
ements of the diagonal matrix are the Schmidt decompo- 
sition coefficients of the bipartition between the semi-infinite 
chains L(-oo, i) and R(i + 1, oo). This implies that Eq. 
O can be rewritten by = £* =1 A a \^)\i//^), where |i/^> 
and \if/%) are the Schmidt bases for the semi-infinite chains 
L(-oo, i) and R{i + 1, oo), respectively. For the biparti- 
tion, then, the von Neumann entropy S can be defined as 67 
S = -Tr[g L \ogg L ] = -Tr[g R logg R ], where g L = Tr R g and 
or = Til Q are the reduced density matrices of the subsystems 
L and R, respectively, with the density matrix g = | X PX*P|. For 
the semi-infinite chains L and R in the iMPS representation, 
the von Neumann entropy S is given by 



s =-|]4iog4. 



(7) 



In Fig. [10] we plot the von Neumann entropy as a function of 
A for^ = 32. In the entropy, there are three singular points 
that consist of two local peaks ( A = 0.366 and A = 1.180, 
respectively ) and one discontinuous point (A = -1.0). In 
fact, the singular points correspond to the transition points 
from the string order parameters and the Neel order param- 
eter. It is shown that the von Neumann entropy captures the 
phase transitions. The discontinuity of the von Neumann en- 
tropy indicates that a discontinuous phase transition occurs 
between the Ferromagnetic phase and the XY phase. The two 
singular peaks show that the XY-Haldane phase transition and 
the Haldane-Neel phase transition belong to a continues phase 
transition. Especially, it should be noted that, in our iMPS 
representation, the von Neumann entropy can detect the BKT 
phase transition between the XY phase and the Haldane phase. 

Central charge and universality class.— For one- 
dimensional quantum spin models, in general, the logarithmic 
scaling of von Neumann entropy was conformed to exhibit 
conformal invariance 68 and the scaling is governed by a 
universal factor, i.e., a central charge of the associated 
conformal field theory. In fact, in the iMPS representation, 
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FIG. 1 1 : (color online) (a) Correlation length £ as a function of the 
truncation dimension ^ at the critical points A c i = and A c3 = 1. 185, 
respectively. The power curve fittings f = rjx" yields rj = 0.330 and 
k = 1.306 at A c2 = 0, and rj = 0.023 and k = 1.955 at A c3 = 1.185, 
respectively, (b) Scaling of the von Neumann entropy S with respect 
to the truncation dimension x at the critical points. From the fitting 
function Six) = a + b^°gX> me fitting constants are given as a = 
0.139 and b = 0.218 at A c2 = 0, and a = 0.343 and b = 0.166 
at A c3 = 1.185. From (i) k = 1.306 and b = 0.218 at A c2 = 0.0, 
and (ii) k = 1.955 and b = 0.166 at A c3 = 1.185 in (a) and (b), the 
central charges are determined as (i) c = 1.001 and (ii) c = 0.509, 
respectively. 



a diverging entanglement at quantum criticality gives simple 
scaling relations for (i) the von Neumann entropy S and (ii) a 
correlation length £ with respect to the truncation dimension 

*as 



^68-70 



e CK i 

S ~ — log^, 



(8a) 
(8b) 



where c is a central charge and k is a so-called finite- 
entanglement scaling exponent. Here, 77 is a constant. By 
using Eqs. d8aT > and i%b\ , then, a central charge can be ob- 
tained numerically at a critical point. 

In the iMPS approach, the correlation length £ can be ob- 
tained from the transfer matrix T defined in Fig. [TJ(c). Ac- 
tually, for a given x, the finite correlation length in the iMPS 
representation can be defined as £(y) = 1/ log Ip£oCy)/a*iCy)N 
where the yUo and p.\ are the largest and the second largest 
eigenvalues of the transfer matrix T, respectively. In Fig. QT| 
we plot (a) the correlation length £ and (b) the von Neumann 
entropy as a function of the truncation dimension^ at the crit- 
ical points A f 2 = 0.0 and A C 3 = 1.185. Here, the truncation 
dimensions are taken as x - 8, 16, 25, 32, 50, and 64. It is 
shown that both the correlation length £ and the von Neumann 
entropy S diverge as the truncation dimension ^ increases. In 
order to obtain the central charges, we use the numerical fit- 
ting functions, i.e., S(x) — a + b\ogx an d £0t) = WC K - F rom 
the numerical fittings of the von Neumann entropies S, the 
fitting constants are given as a = 0.139 and b = 0.218 for 
A, 2 = 0, and a = 0.343 and b = 0.166 for A c3 = 1.185. 
Also, the power-law fittings on the correlation lengthes £ give 
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thermodynamic limit 36,72 . It implies that a comparison be- 
tween quantum many-body states in different phases can sig- 
nal quantum phase transitions regardless of what type of in- 
ternal order exists in the states. Thus, as an alternative way to 
explore quantum phase transitions, the groundstate fidelity has 
been used in the last few years 7 18 1 . In contrast to quantum en- 
tanglement, the fidelity is a measure of similarity between two 
states. An abrupt change of the fidelity can then be expected 
across a critical point (in the thermodynamic limit). Based on 
understanding such a property of the groundstate fidelity near 
critical points, several measures have been suggested such 
as FLS 72 , reduced fidelity 74 , fidelity susceptibility 75 , density- 
functional fidelity 76 , and operator fidelity 77 . However, it is 
known that the fidelity susceptibility cannot detect a BKT type 
phase transition 56 ' 78 . Thus, the fidelity approaches have been 
thought to be a model-dependent indicator for quantum phase 
transitions. In order to show that a BKT type phase transition 
can be captured by the FLS approach, therefore, we discuss 
the FLS in the iMPS representation in this section. 

Once one obtains the groundstate as a function of the 
anisotropic interaction strength A, the groundstate fidelity is 
defined as F(Ai,A 2 ) = |<^(A 2 )|^(Ai))l- Following Ref. M 
we define the groundstate FLS as 



FIG. 12: (color online) Fidelity per lattice site (FLS) d(A\,Az) as a 
function of the anisotropic interactions A! and A 2 for the truncation 
dimension ^ = 32. The pinch points are seen at (a) A = —1.0, (b) 
A = 0.366, and (c) A = 1.180. These pinch points in the FLS are 
consistent with the phase transition points from the order parameters 
in Fig. [6] In (a) for the ferromagnetic-XY phase transition, the FLS 
discontinuity indicates that a first-order phase transition occurs at the 
transition point A = -1.0. In (b) for the XY-Haldane phase tran- 
sition and (c) for the Haldane-Neel phase transition, the continuous 
behaviors of the FLS across the pinch points implies that a contin- 
uous phase transition occurs at the transition points. It should be 
noted that the FLS captures a BKT phase transition between the XY 
and Haldane phases in (b). 



the numerical fitting constants as k — 1.306 and i] = 0.330 
for A c2 = 0, and k = 1.955 and r] = 0.023 for A c3 = 1.185. 
As a result, the central charges are given by c = 1.001 for 
A c2 = and c = 0.509 for A e3 = 1.185. Our central charges 
are very close to the exact values c — 1 and c = 0.5, re- 
spectively. Therefore, the XY-Haldane phase transition at 
Ac2 = belongs to the Heisenberg universality class, while 
the Haldane-Neel phase transition at A c2 = 1.185 belongs to 
the two-dimensional classical Ising universality class. 



VIII. FIDELITY PER LATTICE SITE FOR PHASE 
TRANSITIONS 

For a quantum phase transition, the groundstate of a sys- 
tem undergoes a drastic change in its structure at a critical 
point 7 1-73 . In fact, the groundstates in different phases should 
be orthogonal because the states are distinguishable in the 



, a r lnf(Ai,A 2 ) 
lnd(Ai,A 2 ) = hm , 

L— >oo /j 



where L is the system size. The FLS is well defined in the 
thermodynamic limit even if F(A\ , A 2 ) becomes trivially zero. 
From the fidelity F"(Ai,A 2 ), the FLS has several properties 
as (i) normalization d(A, A) = 1, (ii) symmetry d(A\, A 2 ) = 
d(A 2 , AO, and (iii) range < d(A u A 2 ) < 1. Within the iMPS 
approach, the FLS d(A\, A 2 )22. is given by the largest eigen- 
value fj.0 of the transfer matrix T up to the corrections that 
decay exponentially in the linear system size L. Then, for the 
infinite-size system, d(A\, A 2 ) = pq. 

In Fig. [12] the groundstate FLSs d are displayed as a func- 
tion of the anisotropic interaction parameters (Ai,A 2 ) with the 
truncation dimension ;f = 32. In the FLS surfaces, it is shown 
that there are three pinch points in the spin-1 XXZ model. 
Each pinch points correspond to the transition points identi- 
fied from the order parameters. In Fig. Q~2] (a), the fidelity 
undergoes an abrupt change, which means that the first-order 
phase transition between the Ferromagnetic phase and the XY 
phase occurs at the pinch point. It is consistent with the dis- 
continuous entropy in Fig. [TO] In Fig. [T2](b), the pinch point 
corresponds to the XY-Haldane phase transition point. Con- 
trasted to the fidelity susceptibility, the FLS is able to detect a 
BKT transition successfully, which is consistent with the re- 
sult of our von Neumann entropy. FigureQ~2](c) shows another 
continuous phase transition, i.e., the Haldane-Neel phase tran- 
sition. Hence, it is shown that the FLS approach can be ap- 
plied to characterize quantum phase transitions as a universal 
indicato r 72 ' 82 . 
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IX. CONCLUSIONS AND REMARKS 

We have investigated the string correlations in an infinite- 
size lattice of a spin-1 XXZ chain. In order to obtain a LRSO 
directly rather than an extrapolated string order in a finite- 
size system, the iMPS presentation has been employed and 
the groundstate wavefunction of the infinite lattice system 
has been numerically generated by the iTEBD algorithm. It 
was shown that the x- and y-components of the string cor- 
relations decay exponentially in the Neel phase, while they 
show a unique behavior of two-step decaying to zero within 
a relatively very large lattice distance in the XY phase. That 
is, there is no long-range transverse string order in the XY 
phase and the Neel phase. However, in the Haldane phase, 
the string correlations are saturated to finite values for a rel- 
atively smaller lattice distance, which shows clearly the exis- 
tence of a LRSO. Consistently, the Neel order does not exit in 
both the XY phase and the Haldane phase. This result verifies 
that both the x- and y-components of the LRSO are the order 
parameters characterizing the Haldane phase. The estimated 
critical points agree well with the previous results as A C 2 = 
for the XY-Haldane phase transition and A C 3 = 1.185 for the 
Haldane-Neel phase transition. 

Further, the behaviors of the von Neumann entropy and 
the FLS have been discussed at the phase transition points. 
Both the von Neumann entropy and the FLS capture the cor- 
responding phase transition points including the BKT point, 
which is consistent with the results from the string order pa- 
rameter. Consequently, the von Neumann entropy as well as 



the fidelity approach based on the FLS can be applied to char- 
acterize quantum phase transitions as a universal phase tran- 
sition indicator. Moreover, from a finite-entanglement scaling 
of the von Neumann entropy with respect to the truncation di- 
mension, the central charges are obtained as c ^ 1 at A C 2 = 
and c =i 0.5 at A C 3 = 1 . 1 85, respectively, which shows the XY- 
Haldane phase transition at A C 2 = belongs to the Heisenberg 
universality class while the Haldane-Neel phase transition at 
Ac2 = 1.185 belongs to the two-dimensional classical Ising 
universality class. 

Contrary to other approaches, a feature of the iMPS ap- 
proach is that, just from the iMPS groundstate, its critical be- 
havior can be captured irrespective of whether a system has 
a finite excitation energy gap or not because, in principle, lo- 
cal and non-local order parameters can be calculated directly. 
Furthermore, von Neumann entropy and FLS can be used as 
a universal phase transition indicator for quantum phase tran- 
sition in the iMPS representation. Hence, this iMPS approach 
would be widely applicable for capturing quantum critical 
phenomena in one-dimensional lattice many-body systems. 
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